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SUMMARY 


The designs of two hybrid resistojets, one of a thrust level 
of 667 mN (150 mlb*) for primary propulsion, and one of 44.4 mN 
(10 mlbr ) for auxiliary control for use with either hydrogen or 
ammonia are described. The design techniques used to forecast 
the temperature distributions using mathematical modeling by a 
personal microcomputer are outlined. BASIC language is used 
throughout to give the designer rapid interaction and control in 
debugging using the interpretive mode. A compiler is added to 
reduce the run times of parametric design studies. 

For both designs integrating a compact first stage coiled 
heater, a relatively high voltage device, with a concentric 
tubular gas heater offers direct matching of the 28 V terminal of 
the spacecraft system. For primary propulsion the maximum heater 
wall temperature needed may be kept to less than 55 °C over that 
of the required gas temperature at the throat for hydrogen or am- 
monia by use of the concentric tubular final stage. 

The analysis of the 667 mN primary resistojet, at a gas tem- 
perature of 2500 *K, chamber pressure of 3.16 atmospheres for 
hydrogen or ammonia separately, gives the following performance. 
The specific impulses are 828 and 406 seconds with electric 
powers of 3.09 and 1.98 kW, respectively. The heater efficiency 
on a total power basis is .96 and .94, giving overall total power 
efficiencies of .81 and .66. 

The use of recently developed grain stabilized rhenium with 
refractory metal oxides is employed for all high temperature ele- 
ments to significantly enhance their endurance strength. High 
purity alumina is used for insulator seals in the required areas. 

The small size and increased life demanded of the auxiliary 
thruster compared to the primary thruster require a more robust 
hybrid heater. A bicoil is used to achieve the required resis- 
tance and heat transfer required than would be from the scale 
down in size from the primary thruster. The bicoil mesh is 
formed by two concentric coils of opposite rotation to form a 
sturdy self supporting structure without insulators. The wire 
intersection points are spot welded to provide adequate rigidity. 

The analysis of the 44.4 mN auxiliary resistojet at a gas 
temperature of 2500 °K and a chamber pressure of 3.16 atmospheres 
for hydrogen or ammonia, results in specific impulse of 757 and 
375 seconds with electric powers of 244 and 162 watts, respec- 
tively. The heater efficiency is found to be lower on a total 
power basis than the primary design at .88 and 83. This results 
in overall total power efficiencies of .60 and .47. 
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INTRODUCTION 


Studies of large manned space stations (circa 1969) for NASA 
showed that high performance resistojets have important applica- 
tions to attitude and orbit control (refs. 1 and 2). Both 
specifically chosen high performance gases and biowaste gases 
available as propellants were studied. These studies showed that 
not only lower propellant consumption over the long life 
resulted, but there were other important advantages. These were 
minimizing contamination of the spacecraft environment, limiting 
noise, vibration, etc. In order to deliver the high performance 
required to be competitive in these studies, these resistojets, 
unlike the much earlier units, required: 1) higher temperatures 

than early flight resistojets of the 1960 period; 2) reliable 
operation for long periods of time (years); and 3) high energy 
efficiency, that is, low thermal losses and high conversion ef- 
ficiencies to jet energy resulting in thrust. 

The present report is a formal summary of informal specific 
design studies of resistojets using high performance gases, H 3 or 
NH=,, performed for Lewis Research Center and delivered in the 
1982-1983 period for internal use when studies for a current 
manned space station were again evolving. By informal reporting 
procedures including technical conferences, details of the design 
calculations, drawings, fabrication procedures, material choices, 
etc., were submitted. The basis for the designs presented here 
were* evolutions of the concentric tubular type resistojet 
developed and evaluated in 8,000 hour life tests in support of 
the early NASA space station studies (refs. 3, 4 and 5) in the 

044 Newton class and* subsequently* further developments (refs. 
6 7, 8 and 9) in the .667 Newton size. A new hybrid feature is 

introduced to the thrusters for matching its terminal voltage 
directly to the spacecraft main power voltage, assumed to be 28 
Vdc. 


A program written in the BASIC language was used for mathe 
matically modeling the temperature fields, heat fluxes, electric 
current and voltages, and gas dynamics resulting from the design 
perturbations tried. The program was run on a personal microcom- 
puter of as low as 128K RAM capability. Finally* it takes the 
opportunity to update the materials used. 

The latest 667 mN (150 mlb* ) concentric tubular design, 
(ref. 7) built (ref. 8) and experimentally evaluated (ref. 9) 
operated at a terminal voltage of 14.7 volts and 208 amperes for 
an applied electrical power of 3051 watts. Freedom to indepen- 
dently adjust the terminal voltage so as to match spacecraft 
power would reduce the power matching (weight) problem substan- 
tially. Further, the simultaneous reduction in current would 
reduce the spacecraft power conductor size. However, early at- 
tempts to increase this voltage by reducing the wall thickness 
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of the ohmically heated concentric shells or increasing the num- 
ber of such tubes was reported (ref. 6) to not be feasible for 
such a design. The contribution of the current design is a hybrid 
feature to be described , which is a first stage heat exchanger in 
series, both electrically and gas dynamically, with the con- 
centric tubular design. It consists of a multicoil type ex- 
changer. This type design has the classical high wire over- 
temperature characteristic compared to resulting gas temperature. 
This is not a serious problem in this design in that its tempera- 
ture range is low. It works well within long-life criteria at the 
lower temperatures, yet the coil allows adjustment of the overall 
voltage readily to 28 volts when connected in series with the 
concentric tubular stage. The second stage concentric tubular 
design can be made to have the gas finally approach very closely 
to the maximum wall temperature, a feature favorable to high per- 
formance . 

The problem of designing for 28 volts at the terminals in a 
44 mN thruster is significantly more difficult than in the 667 mN 
size. A different kind of additional stage in addition to the 
first stage coil, a bicoil, is required. 

Analytical and experimental studies of viscous and transi- 
tion flow nozzles are required before such designs can be made 
with great confidence. This was beyond the scope of this effort 
and is an area for future study. However the nozzle efficiencies 
used in the projected performance are attainable, and have been 
realized before (refs. 6 and 11). 

The resistojet is simple in concept. There are basically 
two sections: 1) an electric gas heater, and 2) a gas ac- 
celerator or thrust producer (nozzle). The first is primarily a 
heat transfer problem. The second is supersonic - laminar to 
transition flow field problem involving miniature nozzles. 

In order to address the major thermal impact of the addition 
of the first stage heater, the mathematical model herein 
described was developed. While the notion of thermal modeling to 
determine the temperature field is an old one, the one presented 
here in light of modern personal computers allows the design en- 
gineer to begin to obtain results from his adjustments to the 
design. Reference 10 describes a subsequently developed, more 
generalized approach to the problem by the same authors. It is 
considerably more rapid and general than that reported here. 
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DESIGN OBJECTIVES 


Two resistojets were designed. One is for primary propulsion 
and the other for auxiliary. Both were to serve as laboratory 
models usable for research on hybrid propulsion concepts. This 
included the ability to disassemble the units following a test. 
At the same time, flight weight design philosophy was to be used 
as far as feasible. The resistojets were to be operable on 
either hydrogen or ammonia. Performance goals for both thrusters 
with either propellant are shown in table I. The independent 
design variables include a thrust of 667 mN (150 mlb,*) for the 
primary thruster, and 44 mN (10 mlbr) for the auxiliary. Both 
units are designed for a gas temperature of 2500*K, and a chamber 
pressure of 3 atmospheres. An important goal was that they be 
operable on 28 volts with hydrogen and at least 22 volts on am- 
monia. 


While operation at 28 volts is clearly desirable, the dif- 
ficulty of achieving this in the smaller auxiliary thruster 
designs was recognized. A trade off between thruster design dif- 
ficulties and the alternate requirement for power adaptation 
should be examined. the primary thruster had a minimum design 
life goal of 2000 hours with the capability of withstanding 400 
interruptions of electrical power; while, for the auxiliary 
thruster, a goal of 8000 hours with 4000 cycles was set. 

Numerous considerations enter into the selection of 
materials for a resistojet. Properties of interest include: 
electrical resistivity, ductility after high temperature opera- 
tion, low and high temperature strength, weldability, and for- 
mability to name but a few. Since life requirements are high, 
the design task tends to focus on achieving the required life 
while maintaining reasonable cost and weight figures. 

The concentric tubular type resistojet was first introduced 
in 1962. Because of early experimental work, criteria for es- 
timating the life of this type of unit are fairly well known 
(ref. 3). The materials used for construction, besides serving 
their functions, require resistance to the long term effects of 
high temperatures. Life governing factors are sublimation of 
heater surfaces, redeposition of vaporized material on cooler 
critical areas, endurance stress rupture, creep elongation, 
cyclic fatigue, chemical attack of expellants on walls and in- 
sulators, and metallurgical instabilities. 

An attractive design is one which works each element equally 
based on its limiting criteria. To do this, the designer must 
have accurate data on the endurance of the materials used; and 
have a suitable method for forecasting temperatures, stresses, 
etc., in critical parts. One goal of this program was to 
develop a rapid mathematical computational model to determine 
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effects of thruster design changes on temperature distribution 
heat fluxes and fluid flows in advance of testing. 

The design effort for both thrusters was directed at the 
following areas: 

1 ) Improving and simplifying the space power system to 
resistojet matching by increasing the resistojet's design ter- 
minal voltage to that of the system; 

2) Improving the thruster performance by increasing the 
exit gas total temperature to as close to the heater maximum wall 
temperature limit as possible; 

3) Identification of recommended materials and processes; 

4) Documenting the basis for the design selected by means 
of detailed drawings, materials specifications, fabrication 
processes, and assembly procedures. 


THERMAL ANALYSIS PROCEDURE 


The overall performance evaluation parameters used here are 
defined in Appendix A. The mathematical thermal modeling of both 
thrusters are based upon specific, explicit node diagrams con- 
structed for each case as are shown as figures 1 4 and 1 5 of Ap- 
pendix B. A description of the heat balance program is given. 
Both the general outline of the program and details of the cal- 
culation methods are described there for the two thrusters. The 
reader is referred there for details. 

A later, more generalized, method is described in reference 
10 where the model basis for the thermal analysis is a general- 
ized system of orthogonal elements in cylindrical coordinates. 
This model can quickly be shaped to fit an axisymmetric problem. 
The relationships between elements are activated with references 
set to material resource data banks to use. The initial time 
saved is significantly reduced by this later method. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 

PART I - DESIGN OF A RESISTOJET FOR PRIME PROPULSION 


Description of the Configuration 

The hybrid resistojet concept is described in figure 1 with 
the radial scale exaggerated for clarity. There are three gas 
heating zones in series, the regenerator, the coiled wire heater, 
and the classical concentric tubular heater. The coiled wire 
heater and the concentric tubular heater are separately powered 
for laboratory flexibility. In a flight design these would be in- 
ternally connected. The design contains two vacuum jackets, one 
each between the three heaters, to reduce gaseous thermal conduc- 
tivity losses radially. To prevent high thermal stresses, an ex- 
pansion compensator is designed so as to permit the inner sec- 
tions to axially expand relative to that of the case without 
building severe stresses. Both fibrous and radiation shields are 
used, the locations of electric insulators are noted. 



Figure 1.- Hybrid resistojet concept 
for prime propulsion. 


The resultant geometry of the hybrid 667 mN (150 mlb* ) 
resistojet is shown to scale in figure 2. The callouts of the 
components are identified in table II. This is the final con- 
figuration resulting from the study. The basis for selecting the 
specific geometry is described herein. 
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Figure 2.- Preliminary design assembly of the 667 mN (150 mlbf) 
resistojet for prime propulsion. 


TABLE II 

COMPONENTS OF THE 667 mN (150 mlbr) RESISTOJET 


,( 0# Description 

0. Thermal regenerator 

1 . Hybrid coiled heater-exchanger 

2. Outer pressure case 

3. Inner pressure case 

4. Entrance concentric tubular heater 
pressure vessel (domed end) 

5. Second concentric tubular heater 
(funnel end) 

6. Third concentric tubular heater 
(domed end) 

7. Final concentric tubular heater 

8. Nozzle throat and initial expansion 

9. Nozzle - final expansion 

10. Inner case stem and power terminal 

11. Radial struts - electric conductor 

12. Thermal insulation - low temperature 

13. Thermal insulation — high temperature 

14. Thermal radiation shields 

15. Propellant inlet line 

16. Electric insulators 

17. Coiled heater support 

18. Adapter ring 

1 g ^ Ceramic electric insulator - seal 

20. Expansion compensation bellows 

21 . Thruster mounting pad site 

22. Insulation canister assembly 

23. Electric insulator - bearing 

24. Outer case end 

25. End mount 

* GSRe is grain stabilized rhenium (Ref. 12) 


Material 

CRESS 

♦GSRe/Al^O* 

GSRe 

GSRe 

GSRe 

GSRe 

GSRe 

GSRe 

GSRe 

GSRe 

Mo 

GSRe 
AI3O3 
( fibrous ) 
Zr0 2 

( fibrous ) 

GSRe 

CRESS 

Al 2 0:3 

Alz-O:, 

CRESS 

Al3.0 3 

CRESS 

CRESS 

CRESS 

BN 

GSRe 

CRESS 


Gas h eater . - Physically the hybrid gas heater component is 
divided into three zones, each thermally separated radial y rom 
the others by an effective vacuum Dewar and other insulations. 
In earlier resistojets (ref. 3), with thrust of 44 mN (10 mlbr), 
this technique was shown to be a most effective means of msu a 
tion. A significantly lower gas thermal conductivity is created 
by lowering the Knudsen No. to a number much less than . . 1S 

involves the use of low gas pressure and small characteristic 
lenaths It replaces the continuum regime where thermal conduc 
tivity is independent of gas pressure with the molecular regime 
and is where thermal conductivity is proportional to pressure 
thus controllable. 
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Two such evacuated regions are used. The first is between 
the regenerator (0) and the first powered stage, coiled heater 
exchanger (1). It is filled with fibrous insulation and exposed 
to space vacuum. The second is between the coiled heat exchanger 
and the outer concentric tube exchanger (4). It uses radiation 
shields (14) to reduce the radial radiation term, and is vented 
to space vacuum by means of the inner case stem (10). 

Regenerator: This is the first gas pass. It is located be- 
tween the inner (13) and outer (12) insulation surrounding the 
heater. The gas is heated regeneratively in this pass recovering 
most of the heat that would otherwise have been lost through the 
insulation. Contrary to earlier designs, the gas feed tube en- 
ters at the nozzle end and flows back toward the mount of the 
thruster. This arrangement had been avoided previously to keep 
the higher temperatures away from the mounting flange which tends 
to have a high thermal loss. However, it was found by means of 
the mathematically computationally modeling of the heater that 
the axial distortion of temperature in the heater itself is sub- 
stantially reduced if the first electrically heated pass (in this 
case the second pass, the hybrid coil heater (1)) flows toward 
the nozzle. As a result, the present hybrid design has an es- 
timated gas temperature at the nozzle throat that is only 55 *C 
below the maximum wall temperature. This compares with more than 
110°C in earlier designs, as analyzed here and confirmed ex- 
perimentally in reference 9. The heat loss could be minimized 
by thermally isolating the mounting flange. A further refinement 
that could be included later would involve the use of two 
regenerative passes. The first pass would flow toward the nozzle 
after effectively cooling the housing, while the second pass 
would correspond to the one used in the present geometry but en- 
tering the annulus internally. 

Hybrid coiled heater: This is the second gas pass. It lies 
in an annular passage between two heretofore relatively heavy 
tubes (2 and 3) that served as an earlier regenerative pass and 
pressure vessel (ref. 7). The thickness here has been reduced in 
half from earlier designs. Although these walls conduct current, 
their resistance is low and the electric heating in the walls is 
small. Rhenium wire coils are provided in his pass only. They 
supply most of the gas heating. Alumina guides and liners are 
used to assure that the coils do not touch each other or the 
walls resulting in an electric short. The use of the high tem- 
perature electric insulator, alumina, determines the maximum 
value of the wire temperature that may be used. It is limited to 
2140°K which leads to a design with no coils used after the 
second pass. 
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By sizing the coil wire correctly, it is possible to find a 
design with the heater coil and the nested tubes connected in 
series operating at 28 volts overall. This can be done without 
going to excessively thin wires. In the present design, the wire 
diameter was taken as fixed 1.3 mm ( . 051 inches ) • For the 
laboratory model, this stage was elected to be powered separately 
for diagnostic reasons. Table III summarizes the first stage 
heater. 


TABLE III 

FIRST STAGE HEATER COIL 
PHYSICAL DIMENSIONS 

Heater wire diameter, mm 1.3 


Number of parallel circuits .... 3 

Coils per circuit 2 

Total heater wire length, 

each circuit, mm 880 

No. of turns per coil 5.77 

Coil assembly wire spacing to wire 

diameter ratio 2.5 

Pitch diameter of coils, mm ... 23.5 

Resistance, H 0.1818 


An alternate design for the heater coil involves packing the 
coiled wire tightly in alumina powder which is encapsulated in a 
molybdenum tube. Only the spacing mandrel would be required. No 
electric insulation in the first electrically heated gas passage 
would be required, allowing more flow and room for the encapsula- 
tion of the coil. The latter might be a more satisfactory design 
in that the inevitable diffusion of rhenium over the surfaces of 
the coiled heater support (17) might cause ultimate shorting. 

Concentric tubular heater: Gas passes 3, 4 and 5 are an- 
nular passages between conventional nested tubes. The thin tube 
walls, which are of rhenium, serve as resistance heaters, and 
provide for some additional heating of the gas. This stage is 
described in reference 7. The purpose of the nested tubes is to 
increase the gas temperature above that obtainable with the coil 
heater which is limited by the need for electrical insulators. 
The merit of the concentric design is that the maximum gas tem- 
perature asymptotically approaches closely that of the wall 
producing a high gas temperature. For this reason, it is the 
final stage. The large number of tubes provides the electrical 
resistance required for this stage without using tube walls as 
thin as would otherwise be required. 
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In the sixth and final pass (7), the gas flows through a 
circular tube which connects directly to the nozzle (8). This 
pass is characterized by, a) high thermal conductive losses at 
the nozzle end which cool the wall, b) low convective heat trans- 
fer between the wall and the gas, and c) low resistive heating in 
the relatively thick walled center tube. These characteristics 
are considered fortuitous. The reduced nozzle temperature helps 
to reduce thermal losses, particularly nozzle radiation losses, 
while the low convective heat transfer permits the exit gas tem- 
perature to be nearly as high as the maximum wall value attained. 
The maximum wall temperature naturally tends to be at the center 
of the resistojet as shown later. 

The thruster is thermally insulated by a layer of fibrous 
zirconia insulation (13) surrounded by a layer of fibrous alumina 
insulation (12). Most heat losses through this are recaptured by 
the thermal regenerative assembly. The balance passes outward o 
the (22) insulation canister assembly, hence, thermally to the 
environment through radiation, gas conduction and metal conduc 

tion. 

Note that thermal radiation shielding (14) is shown around 
the nozzle extension. The thermal radiation shield is made in 
two sectors so as to fit around the nozzle extension prior to in- 
“a llation of the insulation and its retaining insulation 

canister assembly (22). 

Nozzle.- No computational gas dynamical analysis of the 
nozzle resulting in its performance efficiency was undertaken. 
She nozzle chosen for use here was that of the earlier 3kW resis- 
toiet (ref. 6). The nozzle efficiency, reported there was 
measured to be 88%. The chamber pressure of 3.16 at ™° s P£ e ^ s 
chosen here as compared to the 8.8 of reference 6 could be 
oected to have somewhat lower performance. The goal values 
chosen for the study can be expected. The nozzle is defined by 

table IV. 


TABLE IV 

NOZZLE DESCRIPTION FOR 667 mN DESIGN 


Configuration 
Half angle, a 
Area ratio 
Throat (geo.), mm 


Conical 

17 . 8 * 

21 2 

1 .339*. 005 


Electric system.- The two exchanger stages for the 
laboratory unithave separate terminals. This is done so as to 
diagnostically evaluate each system separately and be able 
vary the final design power ratios between the first and last 
stage heaters for both hydrogen and ammonia with but one ex 
perimental coil. 
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The electric description of the coiled wire for the baseline 
first stage heater-exchanger chosen is shown in figure 3 and 
Table III The heater is formed as a set of six wire coils each 
with a pitch diameter of 23.5 mm, of 5.77 turns , with an overall 
axial length of 110 mm. The six coils are wound, forming a six 
lead screw (nested coils) about the same axis as the thruster. 
Electrically three parallel circuits are used consisting of two 
coils in series per circuit. Each of the six coils begins and 
terminates at 60* circumferentially from its neighboring coil. 
The two coils in each circuit are connected together at the 
nozzle end and to the respective 28 and 7.7 volt busses at the 
mount end . 

The electric current path for the concentric tubular circuit 
is as follows. From the electric connection which includes a 
thermal dam (ref. 7) at the mount (21), the current flows forward 
through the outer case (2) to the outer case end (24), hence to 
the nozzle (8) and through the inner heating element (7). The 
current leaves the inner heating element, through the struts (11) 
to the next element (6) as is made clear by the figure. It then 
passes through the struts (11) and the second element (5). The 
path continues through its struts and the heater tube (4) and, 
hence, along the inner case and the inner case stem (10) to the 
rear power terminal. 


i 



Circuit 1 


Nozzle 

end 


C i re u i t 2 


Circuit 3 


Figure 3.- Electric schematic of first stage heater. 
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Mechanical system .- The inner assembly is centered in the 
mounting flange by a boron nitride insulator bearing (23) through 
which the stem (10) passes. The inner assembly is attached to 
the outer assembly by a stainless steel bellows (20) which main- 
tains a pressure integrity, but permits differential thermal ex- 
pansion. The bellows is attached to one end of a metallized 
ceramic hermetic seal (19). The other end of this composite part 
is attached to the mount (21) to complete the pressure envelope. 
The basis of this type design is given in reference 3. 

Materials . - From the studies to date, rhenium remains the 
outstanding heating element material. Any further improvement 
would be in the area of stabilizing its grain size to protect the 
geometric shapes of the sophisticated parts and their strength 
during the long high temperature operation required (ref. 12). 

A new electric insulator is proposed and has been used in 
the design. This is high purity alumina. It can be expected to 
be useful in such service to temperatures of at least 2150°K. 


Results and Discussion of the Analysis 

Results of the analysis presented in plots of the tempera- 
tures are shown in Figures 4, 5 and 6. Energy distribution 
diagrams are given in Figures 7a and 7b. 

Temperatures . - One purpose of the analysis is to indicate 
and control locations where the temperature may become excessive. 
It is important that the amount that the maximum wall tempera- 
tures exceeds the exit gas temperature be kept as small as pos- 
sible. Some temperature difference is necessary to drive the 
convective heat transfer. The results show a maximum wall tem- 
perature of 2550 ®K occurred at the entrance to the inner tubular 
heater for the case with a hydrogen gas stagnation temperature of 
2500°K at the throat. Some cooling of the gas occurs because the 
nozzle and the downstream end of the tube are cooled by radiation 
from the nozzle surface. Since the gas is cooled slightly in the 
region approaching the nozzle, it must first be heated to a tem- 
perature higher than the required exit value. Temperatures in 
this region are summarized in table V for the two propellants. 


TABLE V 

CRITICAL TEMPERATURES 
MAXIMUM WALL AND GAS AT THROAT 

Temperature Hydrogen Ammonia 
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Maximum wall 

Maximum gas 

Gas at the throat 


2550 

2529 

2500 


2555 

2523 

2500 


Temperature 
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Figure 6.- Temperature distribution of the 667 mN hybrid 
concentric tubular resistojet operating on ammonia 
from the thermal modeling analysis. 









This amount of over- temperature is considered acceptable. It is 
considerably less than was experienced in earlier designs. Two 
features of the present design keep the over- temperature low: 

1 ) Gas flows toward the nozzle in the first electrically 
heated passage. This passage has the wire coil heater and in- 
creases the gas temperature rapidly to about 2000°K, unlike its 
earlier design (ref. 7). This keeps the temperature at the 
nozzle end of outer case (2) high and markedly reduces cooling of 
the nozzle (ref. 7) by heat conducted through the outer case end 
(24) and across the gas passages. Compare Figure 4 with Figure 5 
to see this temperature distribution difference. 

2) The circular inner tubular heater (7) passage used in 
the last pass connecting to the nozzle has less surface area and 
is not as effective as the other passes as a convective heat 
transfer surface. This allows the tube wall near the nozzle to 
be cooler than the gas without seriously reducing the exit gas 
temperature . 

Although any reduction in gas temperature reduces the 
specific impulse attainable, there is no disadvantage to allowing 
the nozzle to run cooler. In the present design, most of the 
surface of the nozzle is several hundred degrees cooler than the 
gas temperature at the throat. This reduces the heat loss due to 
radiation from the nozzle surface to ambient. If all of the 
nozzle surface was held at 2500°K, the radiation loss given by 
the program would be 382.7 watts compared to an actual 95.61 
watts with hydrogen and 84.07 watts with ammonia. Since 
tradeoffs are involved, a detailed study of nozzle losses is 
needed to assist in the choice of suitable nozzle designs and 
suitable wall temperatures in the nozzle region. 

Losses . - Although thermal radiation from the nozzle was 
found to be the most important loss mechanism in the present 
thruster design (being three or four percent of the input power), 
we are not yet in a position to recommend design changes to al- 
leviate this condition. This is because the radiation loss is 
largely dependent on the nozzle design. If we disregard other 
effects for the moment, this loss can be reduced by: 

1 ) Reducing the nozzle area ratio; 

2) Increasing the inlet gas pressure combined with a reduc- 
tion in nozzle exit area; 

3) Increasing the nozzle divergence angle to reduce the 
radiation "cavity" effect; 

4) Reducing the nozzle wall temperature; 

5) Choosing a surface material with lower emissivity. 
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A trade-off can be expected between radiation loss and other 
forms of nozzle loss such as viscous ones. The design choice 
should depend on a more detailed study. 

An indication of the importance of the cavity effect was ob- 
tained by running the program with the nozzle wall temperature 
set at a uniform 2500°K. At this condition, computational 
analysis gives a radiation loss of 382.7 watts. The emissivity 
for the rhenium wall is .328. The corresponding figure for a 
simple disk of rhenium equal in size to the nozzle exit area is 
217 watts, so the cavity effect (like a hohlraum) corresponds to 
an apparent emissivity that is 76 percent higher. 

Although the present design is not rigorously optimized, it 
illustrates the feasibility of reducing housing losses to some- 
thing in the neighborhood of one percent using good insulation 
and a regenerative cooling passage. At this low level the loss 
effect is within the objectives of this design. The added weight 
is the only compromise in performance. 

It was expected that the losses with ammonia would be 
higher, on a percentage basis, than those with hydrogen. The two 
cases were run with the same throat diameter, the same gas pres- 
sure and the same gas temperature at the nozzle. Under these 
conditions, the electric power input with ammonia is about 65 
percent of the value with hydrogen. If the losses were un- 
changed, they would be more than half again as large on a per- 
centage basis. However, the program showed loss percentages that 
are only slightly higher than those shown for hydrogen. This 
results from the assumption that the gas closely approaches the 
equilibrium condition as it dissociates. On this basis, nearly 
half of the ammonia would dissociate in the regenerative cooling 
passage. The result is more effective cooling, less loss through 
the housing, and some reduction in radiation loss from the 
nozzle. The residence time for ammonia in the regnerative pas- 
sage is estimated to be about 200 milliseconds. A calculation of 
the extent to which equilibrium conditions would be approached 
has not been attempted. 

The computer program did not include an estimate of losses 
in the electrical leads; or, in its final form, losses by conduc- 
tion through the mounting structure. Instead, an idealized lead 
loss has been calculated manually, and the assumption has been 
made that the mounting structure can be combined with the 
electrical leads so that the calculated loss accounts for both 
functions. It has been shown that the lead loss is minimized 
when the leads are sized so that there is no net thermal conduc- 
tion (where they attach to the hot terminal on the engine). See 
reference 3 for the thermal dam design concept. With the rela- 
tively high operating voltage and moderate terminal temperatures 
in the present design, it was found that the minimum lead loss 
could be held to something like one-half of one percent. This 
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suggests that the lead loss can be kept small in a practical 
design. 

Calculated losses are illustrated in diagram form in Figures 
7a and b. The overall thermal efficiency including lead losses 
of gas heater and nozzle is estimated to be 95.9 percent with 
hydrogen and 94.4 percent with ammonia. The calculations were 
carried out and presented with a degree of precision that is 
hardly justified by the accuracy of assumptions used in the cal- 
culation of heat flows. This was done to facilitate the examina- 
tion of the effects of small design changes on performance. 

Gas dynamics .- The Reynolds numbers in both gas heater and 
nozzle indicate the laminar regime everywhere. For the coiled 
heater stage from inlet to outlet, the Reynolds numbers based 
upon flow channel diameter were 185 to 70 and 270 to 86 for 
hydrogen and ammonia, respectively. For the concentric tubular 
heater typical Reynolds numbers were 1300 and 1650 for hydrogen 
and ammonia, respectively. At the nozzle throat the Reyolds num- 
bers for hydrogen and ammonia were 2040 and 2600, respectively. 
The Mach number varied from about .003 at inlet to the first 
stage to .25 at the exit just before the nozzle. 

Power supply requirements .- Although an effort was made to 
size the heating elements to simulate a case with the wire coil 
and the nested tubes operating in series, the program was ac- 
tually run with two elements on separate circuits so that tem- 
peratures at the coil exit and at the nozzle could be controlled 
separately. This was a requirement for the initial laboratory 
engine to provide greater operating flexibility. The heat 
balance shown in table Via gives the following operating condi- 
tions (omitting lead losses). 


TABLE VI 

POWER SUPPLY REQUIREMENTS - 
WIRE COILS AND NESTED TUBES 

a) Initial Results for Independently Powered Stages 


Propellant 

Element 

Voltage 

(V) 

Current 

(A) 

Power 

(W) 

Resistance 

(fl) 

Hydrogen 

Coil 

Concentric 

19.03 

9.59 

102.3 

117.3 

I 946 

I I 25 

.1860 
.081 8 


System 



3070 


Ammonia 

Coil 

Concentric 

14.52 

7.51 

87.9 

92.2 

1276 

693 

.1652 
.081 5 


System 



1969 
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Since the two stages must operate electrically in series, 
their design currents must be the same. The resistance of each 
stage is adjusted uniformly to give the same i^R as before. This 
requires a minor adjustment in wire diameter and tube thickness 
which has no effect on the validity of the thermal analysis. 
TABLE VIb summarizes series operation. 


TABLE VI (Continued) 

b) Adjusted Characteristic for System (Series) Operation 


Propellant 

Element Voltage 

(V) 

Current 

(A) 

Power 

(W) 

Resistance 

(a) 

Hydrogen 

Coil 

Concentric 

17.74 

10.26 

1 09.7 
109.7 

1945 

1125 


.1620 

.0935 


System 

28.00 

1 09.7 

3070 


.2555 

Ammonia 

Coil 

Concentric 

14.26 

7.74 

89.5 

89.5 

1276 

693 


.1593 

.0865 


System 

22.00 

89.5 

1969 


.2458 

Estimated 

performance 

. - The 

estimated 

performance 

presented 

in Part I is summarized in 

table 

VII. 





TABLE VII 

ESTIMATED PERFORMANCE FOR PRIMARY PROPULSION 


Propellant 


h 3 

NH 3 

Thrust, mN 


667 

667 

Gas temperature, 

•K 

2500 

2500 

Chamber pressure, 

atm 

3.16 

3.16 

Specific impulse. 

sec 

828 

406 

Expellant flow, mg sec -1 
Electric power, W 

82.2 

3085 

167.5 
1 978 

Total power, W 


3438 

2078 

Terminal voltage, 

V 

28 

22 

Total power effic 

iencies 


0.626 

Overall, n 0 


0.774 

Heater, r 'n 


0.959 

0.944 

Nozzle, r iis< 


0.807 

0.663 
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PART II - DESIGN OF A RESISTOJET FOR AUXILIARY PROPULSION 


The Problem of Scaling Down 

A resistojet designed for 44.4 mN (10 mlb<- ) of thrust is 
significantly smaller than the 667 mN (150 mlbr ) size considered 
in Part I. If it were feasible to simply scale down the Part I 
design, all dimensions would be reduced by a factor of about 
3.87. However, a number of difficulties would arise if this were 
attempted . 

1 ) The requirement for a long life for the auxiliary 
thrusters 8000 hours versus 2000 for the primary, generally 
places a lower limit on the heating element thickness; which, for 
the auxiliary thruster, is incompatible with the scaled down 
dimensions . 

2) It would be desirable to operate the resistojet directly 
connected to a standard 28 volt power supply. A high operating 
voltage decreases the lead losses as well as avoiding power con- 
ditioning losses. 

It can be shown that, although direct scaling of a heating 
element to smaller size increases the resistance, the increase is 
not enough to maintain the same operating voltage. It is clear 
that a new type of heater coil geometry is needed to provide the 
necessary life, while avoiding the losses and the added weight 
associated with low voltage and the use of power conditioning 
equipment . 

3) As the size of a resistojet is reduced, the thermal 
losses due to conduction become an increased fraction of the to- 
tal power. This is because the temperature gradient in the in- 
sulation (and in other conductive heat flow paths) increases as 
the size is reduced. Even with careful design, it appears to be 
necessary to accept lower efficiencies in the ten millipound size 
range . 


Description of the Configuration 

The key features of the final configuration are shown in 
figure 8 and table VIII. The detailed design features such as 
metal to ceramic seals mounting pads are not shown in figure 8, 
but are similar to those of figure 2, the 667 mN thruster. It is 
significantly smaller than its predecessor which is described in 
reference 3. Figure 9 shows these two units in profile outline 
for comparison purposes. The reduction in size is largely at- 
tributable to the use of coil type heating elements to augment 
the heating action of the concentric tubes. Particularly in small 
size engines, the total tube length becomes excessive if the full 
28 volts is required across the concentric tubular heaters alone. 
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gure 8.- Preliminary layout of hybrid 44.4 mN (10 mlbf) 
resistojet for auxiliary propulsion. 


TABLE VIII 


COMPONENTS OF THE 44.4 mN (10 MLB r ) AUXILIARY RESISTOJET 


No. 

Description 

Material 

0 . 

Thermal regenerator 

CRESS 

1 . 

Entrance coiled heater exchanger 

GSRe/AlarOo. 

2 . 

Bicoil heater-exchanger 

GSRe 

3. 

Inner pressure case 

GSRe 

4. 

Outer pressure case 

GSRe 

5. 

Entrance concentric tubular heater 



(domed end) 

GSRe 

6 . 

Final concentric tubular heater 

GSRe 

7. 

Nozzle throat and initial expansion 

GSRe 

8 . 

Nozzle - final expansion 

GSRe 

9. 

Inner case stem 

Mo 

10 . 

Thermal insulation - low temperature 

Ala-Oa 
( fibrous ) 

11 . 

Thermal insulation - high temperature 

Zr0 2 

( fibrous ) 

12 . 

Thermal radiation shields 

GSRe 

13. 

Propellant inlet line 

CRESS 

14. 

Electric insulators 

AI2O3 

15. 

Coiled heater support 

Al ^0 3 

16. 

Ceramic electric insulator - seal site 

Al^O^ 

17. 

Expansion compensation bellows site 

CRESS 

18. 

Thruster mounting pad site 

CRESS 

19. 

Insulation canister assembly 

CRESS 

20 . 

Electric insulator - bearing 

BN 

21 . 

Outer case end 

GSRe 


Gas heater.- The auxiliary thruster design differs in 
several respects from the prime thruster design previously 
described. These changes are intended to alleviate the problems 
associated with small size engines. 

1 ) Two passes are used for the regenerative preheater in 
an effort to recover more of the heat lost by conduction through 
the thermal insulation. 

2) The heater coil and its annular passage (the third 
pass) are truncated in length at the stem end to reduce the flow 
of heat from the heater to the mounting base and the stem. Con- 
vective heat transfer benefits from the high temperature 
gradients in small engines. The remaining passage length appears 
to be adequate for initial heating of the gas. 
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b) Auxiliary thruster of 
this study 


Figure 9.- Comparative profiles of the auxiliary thruster 
(44.4 mN) of this study with its predecessor. Model II 
(10 mlbf) life test unit (ref. 3). Dimensions in mm. 




3) The fourth pass is an annular passage enclosing a 
tubular mesh heater. The mesh is formed by two concentric coils 
of opposite rotation. The wire intersection points are spot 
welded to provide adequate rigidity while maintaining high 
electrical resistance. The heater structure is rigid enough to 
allow the coil to be free standing without the need for insulat- 
ing liners on the tube surfaces. This pass can therefore run at 
the maximum allowable rhenium temperature. It is a highly effec- 
tive heater. It takes the place of three or four passes of 
nested tube heaters, and does the job with higher electrical 
resistance permitting direct operation from a 28 volt power 
source. The fifth and final pass consists of the center tube 
which carries electrical current and provides a small additional 
amount of heating. The action is similar to that of the final 
pass used in earlier designs. Some heat is lost due to the cool- 
ing effect of the nozzle, but the level of loss is kept moderate 
by designing for reduced heating in this pass. 

4) Simple hand calculations were used to select a modified 
stem design that could be expected to run with a satisfactory 
temperature in the seal region. the selected design increased 
the stem length to three centimeters and added radiation fins at 
the extreme end of the stem. The fins were taken as equivalent 
to a disk of alumina 5 centimeters in diameter and 0.25 cen- 
timeters thick. The presence of the fins increases the thermal 
loss through the stem, but the loss is held to less than 8 watts 
by the increased stem length. The net effect of the two changes 
is to increase the temperature drop in the stem. 

5) To further reduce the temperature in the stem, a radia- 

shield was introduced across the cavity joining the interior 

surfaces of the stem with the hot hemispherical ending of the 
concentric tubular stage. 


Nozzle . - The nozzle is defined by table IX. 


TABLE IX 

NOZZLE DESCRIPTION FOR 44.4 mN DESIGN 


Configuration 
Half angle, a 
Area ratio, (geo. ) 
Throat (geo.), mm 


Conical 

18° 

1 78 

. 348± . 005 


Systems. - The specific characteristics of the heating sys- 
tem are described by figure 8 and table X. The materials of con- 
struction are the same as used for the 667 mN thruster. 
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TABLE X 

PHYSICAL DIMENSIONS 


a) Coiled heaters 

First stage Second stage 

( 1 ) ( 2 ) 


Heater wire diameter, mm 
Number of parallel circuits 
Coils per circuit (series) 
Circuit wire length, mm 
Total heater wire length, mm 


.254 

3 

2 

61 

365 


.203 

4* 

1 

148.6 

595.4 


*2 clockwise/2counter 


b) Concentric tubular heaters 

Domed Pressure 
Vessel (5) 


Tubular 
Inner ( 6 ) 


3.85 1.35 

Outer diameter, mm .125 

Wall thickness, mm * 1 


Results and Discussion of Results 

Estimated performance.- Results of the s ^ te f £° mp J^ th 

analysis of the 44 mN thrus er are a as stagnation tempera - 

hydrogen and ^ “ m °^ a n pP °? e throat in both oases. Temperatures 

are' plotted in figures 10 and 11 ^hele"! igurL 

^^"“hSsff^he^ST 1 ; thruster show a ^^ificant penalty in 

performance associated with -all^engrne^srre. ^ ^ t 

heater efficiencie Derce nt to 82 percent for ammonia. As 

for ^frogen^ and fr P adially outward through the insulation 

expected, tne near uu t-ntal electrical power in the 

is a much larger fraction of total^electr^ s £ ecial care in 

small engine case, maki g . Dass aaes to recover most 

the design of the regenerative heating passages to 

of this energy. 

• another P respect" Ce £ Sxi^ ^nSe^u^" SlS |oSS 

W IBS'-, ra =• =» 
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Figure 10.- Temperature distribution from the thermal modeling analysis 
of the 44.4 mN hybrid concentric tubular resistojet operating on 
hydrogen . 
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Figure 11.- Temperature distribution from the thermal modeling analysis 
of the 44.4 mN hybrid concentric tubular resistojet operating on 







Larger portions of the heat in the gas are lost by convec- 
tion to the cool wall near the nozzle; and in turn, the wall is 
cooled to a greater extent by increased conduction into the hous- 
ing and by increased convection to the gas in the adjoining pas- 
sage. However, the latter effect has been accentuated by the use 
of the second heater coil which heats the gas rapidly resulting 
in a strong axial temperature gradient. Unfortunately, the cool 
end of the pass is adjacent to the nozzle and adds to the cooling 
of the center tube where the flow approaches the nozzle. This 
suggests that adjustments might be considered in the heating ele- 
designs to increase the heating rate near the nozzle. A 
trade-off between high specific impulse and good efficiency is 
involved. 

Table XI summarizes the projected performance of the 44mN 
design resulting from this analysis. 


TABLE XI 

ESTIMATED PERFORMANCE FOR AUXILIARY THRUSTER 


Propellant 

II 3 

nh 3 

Thrust, mN 

44.4 

44.4 

Gas temperature, °K 

2500 

2500 

Chamber pressure, atm 

3.16 

3.16 

Specific impulse, sec 

757 

375 

Expellant flow, mg sec 1 

5.99 

12.10 

Electric power, W 

243.6 

161 .9 

Total power, W 

269.4 

169.1 

Terminal voltage, V 

28 

22 

Current, A 

9.053 

7.351 

Total power efficiencies 
Overall, n 0 

0.597 

0.467 

Heater, n M 

0.884 

0.823 

Nozzle, 

0.673 

0.567 
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CONCLUDING REMARKS 


1 . The design technique used to forecast energy and tem- 
perature distribution in resistojets using mathematical modeling 
by small microcomputer was found to rapidly predict results for 
the designs here. 

2. Where compared with earlier experimental tests the 
model did explain earlier results, resulting in their design cor- 
rection. 

3. The technique of raising appreciably the voltage 
characteristics of concentric tubular resistojets for better 
matching the space station power sources by using a hybrid design 
of a high resistance inlet stage in series was successful. The 
28 volt terminal voltage objective can be met for all designs 
studied here. 

4. The exit gas temperature can be made to closely ap- 
proach maximum wall temperature by use of the modeling technique. 
The differences ranged from 55°C in the 667 mN design. It is 
more difficult with the smaller 44.4 mN thruster where the maxi- 
mum wall temperature required was 186°C over that of the gas at 
nozzle throat. 

5. The study of gas dynamic performance in nozzles by 
mathematical modeling not undertaken here is recommended for fu- 
ture study. 
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APPENDIX A 


DEFINITION OF PERFORMANCE PARAMETERS 

Specific impulse , I»,», seconds 

I — = F/m ( A1 ) 

where F = measured thrust by dynamometer, q r 
ft = propellant mass flow, g-sec -1 

Electric power to terminals , P E , watts 

P e = Et x I* ( A2 ) 

E-b = electric voltage difference between 
thruster terminals, V 
I* = electric current through terminals 

Heater Resistance , R, ohms 

R = E^/Ib, (A3) 

Efficiency - total power overall , n 0 

Mo = F x I,„/(20.8 (P= + Pd) (A4) 

Efficiency - electrical power overall , n 0 * 

Mo* = F x I,„/(20.8 Per) ( A5 ) 

Initial power in gas , P t , watts 

Pi = iii h 4.186 (A6) 

h = absolute enthalpy, ideal gas, cal/gm 
Power in the jet , P_i, watts 

P , = P^ + Pi - Pl_ ( A7 ) 

Pi_ = heat lost from the thruster prior to 
the exhaust jet, W 

Heater efficiency , 'i h 

= P. J /(P= + P d ( A8 ) 

Nozzle efficiency , r ' rn 

n n = F x I„_/ (20.8 Pj) ( A9 ) 
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appendix b 


basis of mathematical computational model 
B y L. Barker 


An IBM P.C. running a to 126 K bvtes o£ memory and 
resistojet. The compiler uses 2 v P t o run. This permitted a 
then requires a total of abo t 192 nodes. Because the total 
node model with about ^ the BASIC chain command 

£ TuTeall combined or three subprograms. 


The Model 

A program has been developed for ^timating ^te*^ types 

temperatures and heat f ^^resented by a network of nodes 
described. The thruster el P ments> sinC e the thruster util- 

joined by heat transm g nassaqes, the node diagram is 

izes a large number concentric There are 22 4 nodes 

relatively complex. S g nea rly 500 heat flow paths. The 

or junction points ejected y ^ re next ass igned, which al- 

material identifiers o mater ial properties subroutine banks as 

lows access to the needed ma r 1 P P wit h temperature, emis- 
reguired, e.g., electrical condi uctiv^y linkages bet ween the 

sivity, etc. The various embodied in the program. Heat 

various adjacent elements are diacent no des and nonasym- 

transfer procedures for tre ^^" g an Individual basis. The net 
metric components are prov deter mined for its assumed tern- 

heat flow into an element can be ge t H eat flows between 

perature and that of its ^acen^^ adjusted suc _ 

nodes are estimated, an satisfactory heat balance is 

cessive approximations unti t he gas temperature is 

achieved at all nodes. At t and the electrical cur- 

estimated at each node in . th ® p f ^ obtain the desired heater coil 
rents are adjusted successively the desired gas 

temperature at the exit to t throat. When the net heat 

stagnation temperature at the noz It brium temperature 

- — “ 


The Method 

To solve for a steady ,-state f temperature ^ribution^the 

detailed procedure calculated using assumed values 

heat flows between ? 1 J h "° d ^ e a aumme a for each node. If the sums 

S not^suf f iciently SSirShT’Sper.ture. are adjusted and the 
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calculations repeated. The 
used is a fraction (usually 


systematic temperature correction 
.7) of the following expression: 


N 

- 2 q (J 
j=1 

AT, = 

N 3 q , j 

Z 

j=1 3T, 

where q, j = the heat flow into node i from node j. 

There is a test for convergence. This section requires that 
heat balance error be no more than 0.1 watt for any node and 
that the heating coil temperature at the down stream end and the 
gas exit^ temperature differ from the prescribed values by no more 
than 0.1 *K. If the conditions for convergence are met, the pro- 
cedure leaves the main loop and jumps to the printout section, 
and if not, the iterative process continues. 


next approximation equation for temperatures above, 
it has been found that the fraction (.7 used here) must be less 
than unity to stabilize the convergence process. Too high a 
value can result in increasing errors of alternate sign, while 
too low a value will reduce the rate of convergence. 


When the temperatures have 
and temperatures are printed out 
on a dot matrix printer. For a 
proach is used, except instead 
tion, a similar expression i 
specific heat data for each node 
time interval. The models used 
follow. 


converged, values of heat flows 
and the temperatures are plotted 
transient solution, the above ap- 
of the above temperature correc- 
s used which includes mass and 
as well as a maximum allowable 
for various heat flow mechanisms 


Conduction 

The equation used is the Fourier equation: 
3 T 

q = -kA 

3 x 


In this equation heat flow towards 
the program heat flow into a 
direction . 


increasing x is positive. In 
node is positive regardless of 


For a given qeometry one 
that A is a function of x. 


can assume that q is a constant and 
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/ INITIAL CONDITIONS 

/ Introduce dimension statements, 

Define geometry and operating conditions, 
Assign values to physical constants, & 
Read in initial guesses for temperatures, 
and electrical currents. / 


ELECTRICAL CURRENT A DJUSTMENTS 
Adjust current in the heater 
to reduce errors in the heater temperature 
and the exit gas temperature. 

' * , 

CONVECTIVE HEAT TRANSFER 
Estimate convective heat flows and their 
derivatives. Find gas temperatures starting 
at the inlet (with a known temperature) and moving 
downstream one node at a time. (Wall temperatures 
are held constant du ring this calcu lation). | 

* zz 

OTHER HEAT FLOWS 

Calculate resistive heating and radiative 
and conductive heat flows with their derivatives. 
(These may be found in any order). 

Wall temperatures are again held constant. 

ERROR CALCULATION 

Determine heat balance errors at each wall node. 
Find errors for the maximum heater temperature ana 
for the gas temperature at the nozzle throat. 


^XAre errors \ 
within acceptable; 
V limits? / 


NEXT APPROXIMATION 
Adjust wall node temperatures to 
reduce heat balance errors making 
use of the error magnitudes and 
the heat flow derivatives. 


PRINTOUT 


Figure 13.- Flow diagram for the 
Computational Mathematical Model 
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Figure 15.- Node diagram - heater 
computational mathematical model 
44.4 mN auxiliary thruster. 













Or, rearranging terms: 


q dT 

= -k(T) 

A(x) dx 

For axial conduction A(x) = A, a constant. 


and q 


A 


(x 2 - Xx ) 



( T)dT 


For radial conduction A(r) = 2nr dx 


and q 


where , 


2ndx 



I* o — r o >_-i m* r 

r i — r 1 nn*r 


T i 

j k ( T ) dT 


For nodes that are neither axial nor radial, a linear area 
variation is assumed approximately orthogonal to flow path. 


Thus , 


and q 


A ( L ) = aL 


a 


In Lj,/L i 


Tx 


/ 


T* 


k(T)dt 


To simplify the problem, when adjacent nodes are not of the 
same material, a heat flow resistance approach is used. 

2 ( Ti- - Tx) 

In the program, qx = 

(Rx + R 2 ) 

Where = dx / ( Ak ( T ) ) 

R r .d = ln(r 0 /rx) / (2ndxk(T)) 
and R«»n = lnfLz/Lx) / ( ak ( T ) ) 
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Convection 


The equation used for tubes is: 
q = (ANuk / D h ) (T w - T„) 

Where D,-. = hydraulic diameter = 4(flow area/flow perimeter) 

If the propellant temperature is not known, then the energy 
received by it is equated to its heat capacity to obtain 

T„ = Tw + (T„o - T„)e < - hn ^ c ’” 

where T„ 0 = initial gas temperature. 

and then 

q w - (T„ 0 - T„ ) mCp 

The approach for annual passages is described in Kays (Ref. 29). 


Ai 

Where q i = (K*Nui i ( T i -T m )+8 iK 0 NUoa(To“Tm) ) 

A 0 

q 0 = — (0 o k,Nui i ( T i -T m ) + k 0 Nu a o ( T 0 T m ) ) 

D(~, (1- 0i©o) 

In the nozzle high velocity region, the approach used is 
described in Eckert and Drake. The equation used is: 

q = (ANuk / Dn) (T w - T« w ) 

where k and Nu are determined for the reference temperature 
T* = T_ + . 5 ( T w - T_) + .22(T.w-T-) 

Y - 1 

where T 0 = T_ ( 1 + M-.^ 5 ) = stagnation temperature 

2 

Y -1 

T„w = Tw ( 1 +r M. 2 ) = adiabatic wall temperature 

2 

and the recovery factor r = Pr 1 ' 3 for laminar flow and 
r = Pr 1 ' 3 for turbulent flow. M_ = the mach number input to the 
program as a constant for each node in the nozzle which is calcu- 
lated as a function of gamma (y) and the ratio of cross-sectional 
area to throat area. 
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Radiation 


Given two surfaces, the equation used is 


A i ( et.i- - ee.1 ) 



where e b z = oTj** 

Usually Fi 2 = 1 and this equation reduces to Christianson's 
equation. For radiation to ambient A z -> ® and the equation be- 
comes 

q x = Ax e x (e b2 -ebi ) where e b2 = oT“„ mfc . 

For more than two surfaces simultaneous equations are used 
to solve for the heat transfer. 

For each node, adaptinq the terminoloqy of Eckert and Drake 
(Ref. 34). 

Q = q/A heat transfer per area per time 

et, = oTx. the emissive power, net radiation per 

area per time 

H = net heat flux arriving at node 

B = net heat flux leaving the node 


r 

m ~~ 

(B/e*) 

the apparent emissivity 

£ II = 


the absorbed part of H 

Then, 

Qi = e.Hi 

— ^ 1 0 1 


1 

N 

where 

Hi = 

S ( A^F k 1 Sbk ) 


A, k=1 
N 

or Hi = S (F( k e ««k: ebk) 

k=1 


Thus, for each node i 

Q i = - £ x 


N 

e b i + 2 (Fixe c .k etii< ) 

k=1 
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Also, 


Bi = Hi - Qi . From this, one can derive 
(1 - fi ) N 

£ (ft 1 — £ 1 — £ (Fi k £ « k*; C fc* k< ) 

6b 1 ]c = 1 

This set of equations can be solved for the £ „ ' s and by sub- 
stituting into the equation q, can be found. Or, when £i<> 1, 
the equation can be reduced to 

e i 

Q i — Cbt ( s «, i — 1 ) 

1 -£ i 
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